function [etilde_ij] = entrants(q_ij,s_i,cobb_par,M,delta,sf)
%UNTITLED3 Summary of this function goes here
%   Detailed explanation goes here

G_ij = q_ij./nansum(q_ij,2);
q_i = nansum(q_ij,2);
e_i = inverse_m(q_i,s_i,cobb_par,M,sf);
etilde_i = e_i-delta*(e_i-q_i);
etilde_ij = etilde_i.*G_ij;

end

